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ON  THE  SOLUTION  OF  A  CLASS  OF  DUAL  INTEGRAL  EQUATIONS 

1  .   Introduction 

The  purpose  of  this  paper  is  to  discuss  the  solution  of  a  class  of  dual 
integral  equations  which  appear  in  the  formulation  of  electristatic  and  elec- 
tromagnetic boundary  value  problems  possessing  circular  symmetry.   Two  classes 
of  equations  are  discussed  out  of  which  the  first  one  admits  a  closed  form 
solution.   A  Fredholm's  equation  of  the  second  kind  is  derived  for  the  second 
class  and  iterative  means  of  solution  are  suggested.   Some  particular  advantages 
of  the  formulation  are  pointed  out, 

2.   The  Dual  Integral  Equations  -  A. 

The  equations  under  discussion  are: 


a 

1 


\a  f(\)  J  (Xp)  d\  =  I   H(X)  J  (Xp)  d\  0<p<l       (1) 


f(\)  J  (\p)d\  =0  p  >1  (2) 


where  H(X.)  is  a  given  function  a  >-l-2u,  and  p.^.0. 

1       2 
Equations  similar  to  these  have  been  discussed  by  Tranter   Noble   and 

others.   A  method  for  obtaining  a  solution  of  the  above  equations  in  a  close 

form  will  be  given  in  the  following  section.   The  method  outlined  will  be 

developed  along  a  line  which  initially  is  somewhat  similar  to  that  of  Tranter's 

The  chief  advantages  of  the  present  method  over  that  of  Tranter's  are  the 

following; 

(i)   The  solution  is  obtained  in  a  closed  form  rather  than  in  a  series 
form  given  by  Tranter, 

(ii)   The  present  method  is  especially  suited  for  obtaining  a  Fredholm's 
equation  of  second  kind  for  the  dual  equations  -  B?  whereas  Tranter  's  method 
involves  the  solution  of  an  infinite  set  of  equations  for  this  case. 


The   equations   B   ares 

oo  OO 

I         \a     fl+T(\)]       f(\)    J    (ty»d\  =  |       H(\)    J    (Xp)d\  0<p<l         (3) 

o  ■'o 


f(\)    J    (ty))d\  =0  P  >1  (4) 

o 


Where  T(\)  is  a  known  function  and  T(\)  ->  0,  or  to  a  constant  for  large 
positive  values  of  X.   The  equations  (3)  and  (4)  are  of  considerable  interest 
and  will  be  the  subject  of  discussion  of  a  later  section. 

3.   Method  of  Solution  -  Set  A 

As  a  first  step  assume  that  f(X.)  is  represented  by  the  series 

oo 

f(\)  =  \1_k   V  C   J0     ,,..,      k  >0  (5) 

/_±      m   2m+^+k(\)' 

m=o 
Such  a  choice  suggested  by  Tranter  is  in  view  of  the  relation 

OO 

r 

xlh  J   o   JX>    J  <V»  =  0  P>1        (6) 

v+2m+h       u 


The  representation  of  f(\)  in  (5)  automatically  makes  it  satisfy  (2) .   We 
still  have  to  satisfy  (1) . 

The  first  step  toward  this  consists  of  multiplying  both  sides  of  (1)  by 

p^+1  (i-p2)k_1   ^p  u-ni,  ji+i,  p2) 

and  integrating  w.r.t.  p  in  the  range  0  to  1  where 
(a,Y,x) 


%  =  2F1  ('m>  a+m'  *»> 


=  Jacobi  polynomial  of  order  'p'  and  F 
=  Hypergeometric  function. 


Upon  carrying  out  the  integration,  one  obtains  by  using  the  relation 
(see  Tanter  p.  113) . 


1 

-k 
P  "(J,  x  "r"'  ~^  "u,+2m+kx "  (7) 


p^+1  (i-p2)k_i     0-  J  (\p)  dp  =  \" 


in  (1),  the  following  equation  for  f(\.): 


0  ^0 


Substituting   for   f(\)    from   (5)    there   is    obtained 


[va+l-2kScj  i(\)d\  = 

m   2m+jj.+k         J         2p+jj,+k 


Now  choose  k  such  that 


H(M  X  k  JQ     ,  (M  d\    (9) 
2p+u,+k 


a  +  1  -  2k  =  -1         or       k=l+  ~- ,  (10) 


Since  it  has  been  assumed  that  a  > -1  -2jjl,  it  follows  that 


k>-  -  u.      or      p.  +  k>-  .  (11) 


When  k  is  chosen  according  to  (10),  (9)  becomes 


[X  -1  2  C  j      (\)]  j      (\)  d\  =      H(\)  \  (1+a/2)Jo    ,(X)  d\. 

m=0   m   2m+p,+k    '   2p+jj,+k  I  2p+u,+k 

0  ^0  (12) 

Now  use  the  orthogonality  relation  for  the  Bessel  functions,  viz0, 

X   J„    ,  (X)  J     ,  (\)  d\ 


2m+u.+k      2p+|j.+k  (4m+2k) 


where 


6U  =  o 


for 


b  4   m 


for 


b  =  m. 


and  obtain  from  (12) 


C      = 
P 


H(\)    X~il+a/2)   J  .     (\)    d\ 

2p+|JL+k 


Letting   v  =  jx+k+1    and 


(13) 


S    (\)      =     2       C      «T  ,     (\)    d\ 

v  p=0        p      2p+|jL+k 


we  write   f(X-)    given  by    (5)    as 

f(\)   =  \^~V     S   (X) 


(14) 


Following  Wilkins  ,  a  closed  form  expression  for  S  (. X.)  is  given  by 


Sy(\)  =  H(\)\  a/2   -  \      r  J  t\r)dr       H(t)t  a/2  Jy(tr)  dt    (15) 


provided  the  conditions  specified  by  Wilkins  on  H(\)  \     are  met  and 

v  >  -  — .   The  condition  on  v   is  already  satisfied  because  of  the  condition 

expressed  in  (11).   This  is  also  the  reason  for  imposing  the  condition  on  a 
which  was  stated  in  connection  with  the  equations  (1)  and  (2). 

From  (14)  and  (15)  a  close  form  representation  of  the  solution  is  then 


oo  a 

:(\)  =  x<H--v-a/2>  R{^    _    ^-v  f  p  ^(Xr)dr   f 

1  ^0 


H(t)t"a/2  J  (tr)  dt    (16) 


Recalling  that 

v  =  jjL+k-1   =  |x+a/2 

one  can  rewrite  (16)  as 


f(\)  =  \"U  H(\)  -  X1  Q/2   /    r  J  +a/2(\r)dr       H(t)t  a/2  J  +q/2  (tr)   (17a) 

J1  J  0 

Since  jjl+cl/2  =    v  >-    1/2   one  may   derive   using   Hankel1  s    formula 


H(t)t    a/2  J  /0(tr)dt 

jjL+a/2  J  |jl+q/2 

0 


1  ^OO 

=   H(K)    X-1-72   -        f      r  J   ^(aridr       f       K(tW«/2  J        /2(tr)dt    . 
J0  J0 


Using  this  in  (17a)  one  obtains  an  alternate  form  of  the  solution  as 
*(M  =  [  r  J|x+Q/2(ar)dr   J   H(t)t-Q/2  j^  (tr)  dt         (17b) 

and  one  may  use  one  or  the  other  form  whichever  is  convenient, 

4.   The  Dual  Equations  -  B 

The  equations  under  consideration  in  the  present  section  are  the  same 
as  (3)  and  (4)  and  are  repeated  here  for  convenience.   They  are;   , 


\a    |l+T(\)|  f(\)  J  (\p)d\  = 


H(\)  J  (Xp)  d\    ,         0^p<l    (18) 
0 


f(\)  J  (\p)  d\  =   0  P>1.      (19) 

0 


Following  through  the  same  procedure  which  was  used  to  derive  (9),  one 
now  gets,  assuming  the  following  series  form  for  f(\),  viz., 

f(\)  =  \1_k   2  C  J0    .  (\), 

m=0   m   2m+jjL+k     ' 

the  equation; 


a+1   2k     |1+T(X)|   ?    c 

t  J  m=0      m      2m+jj.+k  2p+u>k  2p+(a+k 

0  ^0 

(20) 

Again  choose  k  =  l+a/2  and  obtain  from  (20) 


[\l      (l+T(\)}  SC   J      ,(\)]  J      .  (\)d\  =     H(\)\  k  JQ     .(\)d\ 
I  L      J  m  m   2m+jjL+k       2p+(j.+k        I  2p+u>k 

0  0  (21) 

Since  (21)  is  true  for  all  p  =  0,  1,  2,  ...  etc.,  one  can  write 

\_1   |l  +  T(\)|     Sv(\)  =  H(\)\"k  +  0(\)  (22) 

where 

S  (\)  =SC  J0    ,(M=2c   Jo   ,   (\),     v  =  n+k-1. 
v  mm  2m+(j.+k      m  m   2m+l+v    '        r 

(j)   (v)  is  a  function  which  satisfies 

r 

0(M  Jo   -,   (Md\  =   0   for  all  p  =  0,  1,  2,  (23) 

/     T  2p+l+v 

Jo 

A  suitable  representation  for  <p   (X)  which  is  suggested  by  the  form  of  (15),  is 

oo 
0(\)  =      J   (\x)  ¥  (x)  dx  (24) 

which  satisfies  (23)  as  is  easily  seen  by  the  substitution  of  (24)  in  (23)  and 
the  employment  of  the  relation 


J        (\)  J^(\x)  dx  =  0   for  x  >  1 


(25) 


It  has  been  assumed  in  the  above  argument  that  it  is  possible  to  change 
the  order  of  the  integrations  involved.   Putting  (24)  in  (22)  we  get  after 
rearranging 


VKK)    ~       1+T(\)        1+T(\) 


J  (\x)  ^  (x)dx 


(26) 


Now,  in  view  of  (25), 


S  (\)  J  (\t) 

v     v 


2C  J   .  „(t)   J  (t)  =  0 

m   2m+l+2       v 


t>l 


Hence,  we  can  derive  the  following  equation  from  (26) 


H(\)\ 


1-k 


1+T(\)     v 


J  (\t)d\ 


1+T(\)   v 


J  (\t)d\ 


J  (\x)  *(x)dx  ,   t  >1.   (27) 


Using  Hankel's  formula  one  can  arrive  at  the  equation 


*(t) 

t 


\J  (\t) 


J  (\x)  ¥  (x)  dx 


t  >1 


(28) 


md  then  rewrite  (27)  as 

OO 


H(\)      x-Q/2  .  ,x    _  J(t)    _  \T(\) 

1+T(\)  V^/2  "      t  1+T<\)    "J140/2 

0  ^0 


J        ,„(U)d\ 


X  J        /2(^x)   ^x)    dx,  t  >  1 

1 


(29) 


8 
Equation  (29)  is  a  FredholmJ s  equation  of  second  kind  for  the  unknown  ^. 
The  equation  of  the  second  kind,  as  is  well  known,  has  the  advantage  that 
iterative  methods  can  be  employed  to  obtain  successive  approximations  of  the 
unknown  and  is  therefore  preferable  to  the  integral  equation  of  the  first  kind 
in  most  instances. 

Some  of  the  particular  advantages  of  the  Equation  (29)  will  be  discussed 
later  in  connection  with  the  examples. 

5.   Examples 

In  this  section  we  shall  illustrate  the  techniques  of  dealing  with  the 
dual  integral  equations  with  some  examples. 

(a)   The  electrified  disk: 

The  first  set  of  equations  to  be  considered  are  the  following: 

f"  f 

x  1  fa)  jo  (xp)dx  =  (2/ 77)  v  I  x~1/2  j      (x)  jq  ap)d\  =  v    for   o<p<  1 

o  Jo 

(30) 

x> 

f(X)  Jq  (AP)  =0  for  p>l 

(31) 

These  equations  are  associated  with  the  problem  of  finding  the  potential 

4 
:  an  electrified  disk  at  a  potential  V  (for  details  see  Tranter  ,  p.  50) . 

The  equations  are  of  class  -  A  and  we  identify  by  comparison  with  (1)  and 


(2), 

a  =  -1,  |1  a  0  and  H(X)  =  (2/ir)  V  J    (A)  X~±/2 

From    (17a)    the   solution   is   given   by 
f(X)    =    (2/TT)    V   J1/2    CmV2    -(2/TT)    V   XV2 


from  Sonines  integral  formulas  (see  Magnus  )  we  have 


f 

DO 

r 

J       r   J-l/2    (Xr)dr 
1 

J    Jl/2(t>    J-l/2« 

o 

(32) 

J1/2  (t)  J_1/2(rt)dt  =0  for  r  >  1. 


9 
The  second  term  in  the  r.h.s.  of  (32)  is  therefore  zero  and  we  thus  have 
the  solution  in  the  simple  form 

few  =  (2/ir)  v  j      an372  (33) 

which  is  in  agreement  with  the  result  given  by  Tranter.   The  use  of  (17b) 
will  of  course  produce  the  same  solution. 
(b)   Coaxial  conducting  disks: 

The  dual  equations  associated  with  thel problem  of  two  charged  coaxial 

conducting  disks  with  radius  a  ^.nd  spaced  apart  at  a  distance  of  2d  have  been 

6  12 

given  by  Love   and  they  have  been  discussed  by  Tranter   and  Noble  .   The 

equations  are: 

x~1a+&~2Xd)fa)jjxp)d\  J~)  v      hgdj  (xp)dx  =  vcp) 

O  I  IT'      I  O 

o  o 


f(A)  j    ap)  -  o 

o 


for 

0<p<a 

(34) 

for 

p>a 

(35) 

where  H(X)  and  V(P)  are  known. 

Comparing  with  (18)  and  (19)  we  identify 


0-    =  -1,    TOO  =  e  2Xd>      n   =;  0   and   H(X)  =  (-)   V 


x  1/2 

Because  (35)  is  to  be  valid  for  P  >  a  ratlier  than  P  >  1,  which  is  true  for 
(19)^  slight  modifications  are  necessary  before  the  formulas  in  Section  4 
can  be  used. 

For  one  thing  the  series  for  f (X)  has  to  be  of  the  form 

f  (X)  -  I1'*       Cm  J2m  ^   +  k  (rt)  =  X1"11  Sv  (Xa)       (36) 


in  order  for  it  to  satisfy  (35) . 


10 


An  equation  corresponding  to  (26)  may  be  derived  and  it  reads 


VaX)  =  1^-   ^1_k   +  17x7*7  Jv  <Xx>    *<x>    dx 


An  integral  equation  for  ^(x)  is  obtained  by  multiplying  both  sides  of 
(37)  by  Jv(^t)  and  integrating  w„r„t.  t  in  the  range  0  to  oo   ,   This  yields 


i+T(^y        dv  v  X) 


X       (Xt)d^ 
1+T(M   ^V 


J  Ax)   *(x)  dx 


which  after  rearranging  may  be  put  in  the  form 


H(X)   x-  2 
1+T('M 


j    (Xt)dX 


*<t) 


'o  K+2 


t>a, 

(38) 


J_1/2(M*, 


t  >  s 


(39) 


Substituting  for  a,    T(^)  etc,  in  (39),  for  the  problem  under  consideration 
we  obtain  the  equation 


(2/77)  V     JWQ  (\a)  J  wo  (\t)  (l+e"2  d)    Xd\ 


1/2 


1/2 


*(t) 

t     "1   1 


^2Xd-   J-l/2(Xt)dX     J-V2(Xx)*(x)dx 


t  > 


(40) 


Equation  (40)  is  the  desired  Fredholm's  integral  equation  for  ^(x) ,   Once  ^(x) 
is  known  for  x  >  at  f(X)  is  found  through  the  use  of  (36)  and  (37), 
We  shall  now  discuss  several  special  cases,  e„g., 

(i)   d  ->  0      (ii)   d  -***>  (iii)   a ->  0,       (iv)   a-»oo  ,  in  connection 

with  (39)  and  (40) 0 

(i)   d-»  0: 

If  d-»  0,  (40)  reads  in  the  limit, 


-<!)v 


Jv2<Va>J-v2(Xt)dX  =  *Jr--5 


X  J.1/2(^t)dX 


J   .  ,„(\x)®(x>dx,  * 

—  1/^ 


(41) 


11 

Since  the  integral  in  the  l.h.s.  of  (41)  is  zero  for  t  >  a,  the  solution 

^(t)  =  o      for  t  >  a 
is  seen  to  satisfy  (41).   From  (36)  and  (37)  we  have 

f  (M  =  ^~k  Sv(a\)  =  \   H(M  \2~2k  =  i  V  J1/2  (aX)  \1/2  (42) 

which  may  be  verified  to  be  in  agreement  with  the  solution  of  (30)  by  compar- 
ing the  equations  (30)  and  (34)  and  their  solutions. 

(ii)   dH>oo 

This  is  a  trivial  case  since  (34)  becomes  identical  with  (30)  for  d  -» oo , 

(iii)   a->  0 

When   a  ->  0,    a   look  at    (34)    and    (35)    will    tell   us   that   f  (X)    should   also   tend 
I   to   zero  for   any   general   H(X) „ 

From    (39),    letting   a  =  0  we   get   after   rearranging   the   r„h,s. 

H(M    ^1/2   Jy    (M:)     (1+T(M~1    d^ 

to  ao  (4.3) 


I    i^TOT    Jv  (Xt)dx 


Jv(Xx)    *  (x)    dx,         t  >  o, 


Using  Hankel  s  inversion  formula  we  get 

-H(\)  X1/2  (1+T(M)_1   =  X(i+T(\))~1    I  Jv(Xx)  *  (x)  dx  (4  4) 

o 
0)  .     From  (37)  one  has  for  a  =  0 

)  Sv(x)  =  0      implying  f  (M  =  0o 

For  small  a,  the  equation  (44)  maybe  used  for  a  starting  approximation  of  ^(X) 
for  the  purposes  of  iteration, 

11  |(iv)   a  •>  »  : 

Again  considering  an  arbitrary  H(X)  and  a  very  large  radius  of  the  disk,  we 
see  from  (37)  that 

»*<*>*  A  xl"h  (4-5) 

Now  since  the  range  of  validity  of  (34)  is  now  the  entire  positive  range 
ty  of  p,  it  is  easily  seen  from  that  equation,  that 

,   |  UNIVERSITY  OF  ILLINOIS 
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l1m)=SvX-h-.  ^g^       ors^f^ 


1+T(X)      ~v  '■  i+T(A) 


for  a  =  x  -1.   Since  k  =  1  +  — ,    (45)  and  (46)  are  seen  to  be  in  agreement. 
Again  (46)  may  be  used  for  a  starting  approximation  for  iteration  when 
a   is  large. 

The  Fredholm;s  equation  (39)  is  therefore  useful  from  the  point  of  view 
of  iterating  both  for  large  and  small  radii  of  the  disks.   This  is  found  to  be 
of  great  advantage  when  applying  the  technique  outlined  here  to  the  problem 
of  electromagnetic  diffraction  due  to  a  circular  disk  or  a  hole  for  both  large 
and  small  diameters.   In  addition  to  this,  through  the  use  of  techniques 
developed  here,  a  Fredholm's  integral  equation  may  also  be  derived  which  is 
particularly  suitable  for  small  radius  a   because  the  kernel  of  this  parallel 
equation  is  small  when  a   is  small.   The  detailed  discussion  on  some  of  these 
topics  will  appear  in  a  later  report. 
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